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Abstract 

A simple and general prescription for evaluating unambiguously the sign 
of the grand-canonical trace of quasi-particle statistical density operators 
(the so-called sign ambiguity in taking the square root of determinants) is 
given. Sign ambiguities of this kind appear in the evaluation of the grand- 
canonical partition function projected to good quantum numbers (angular 
momentum, parity and particle number) in the Hartree-Fock-Bogoliubov 
approximation at finite temperature, since traces are usually expressed as 
the square root of determinants. A comparison is made with the numerical 
continuity method. 
Pacs numbers: 21.60.Jz 
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1 Introduction. 



Nuclei at finite temperature are usually studied microscopically using the finite 
temperature Hartree-Fock (HF) or Hartree-Fock Bogoliubov (HFB) approxima- 
tions (ref. [1],[2]), at least at first level of approximation. The basic quantity is 
the grand canonical partition function or equivalently the grand potential. Such a 
quantity contains an average of all possible contributions from different conserved 
quantum numbers such as angular momentum, parity and particle number. It is of 
course of great physical interest to study the partition function and related ther- 
mal averages using the partition function projected to the exact quantum numbers. 
This is especially true if the projection to good angular momentum is carried out 
without any assumption about axial symmetry, that is, if the full three dimensional 
angular momentum projector is used. As the temperature increases, we expect tri- 
axial shapes to play a role and it is interesting to see how shape transitions are 
obtained at different excitation energies. 

Sometimes the full projector is replaced by a partial projector to good z- 
component of the angular momentum J^. The partition function for a specified 
value of the angular momentum J is then obtained by subtraction between the 
partition function at a Jz and the partition function at J2 + 1. This recipe, however, 
has a basic limitation that requires the exact evaluation of the partition function. 
This limitation is of course a problem when using even accurate approximations 
such as the HFB. Moreover, in the limit of temperature, the HFB ground state 
for a specified value of Jz is not as accurate as the ground state obtained with a 
specified value of J. Also, the above recipe would pose severe problems for odd 
and odd-odd nuclei. Therefore, the use of the exact angular momentum projector 
is highly desirable. 
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In the case of the the temperature dependent HFB approximation, a standard 
result for the trace of the statistical density operator (the exponential of a quadratic 
form in the quasi-particle operators) states that it can be recast as a square root of 
a determinant (ref.[3]). This is a problem if the projected partition function is 
required, since an improper sign can lead to erroneous results. In the past the 
cure for this problem has been given with the so-called continuity argument. This 
argument states that the proper sign can be determined by constructing the statis- 
tical density operator from unity and then by determining the appropriate sign by 
imposing the continuity of the phase of the trace as we progressively rebuild the 
statistical density operator. This has been the recipe followed in ref. [3]. 

Recently this problem has been considered anew using the Grassmann algebra 
(ref. [4]) for the determination of the sign for both overlaps of HFB wave func- 
tions and the trace of the statistical density operator (ref. [5]). Although the results 
obtained were not previously reported in the literature, in the case of the trace of 
the statistical density operator, the sign ambiguity was not fully resolved since the 
vacuum contribution was still left as the square root of a determinant. 

The purpose of this work is to show how all possible ambiguities can be re- 
solved without referring to a numerical continuity argument, which may not be 
easy to implement. In the next section we shall derive the construction of the 
proper sign for the trace of the statistical density operator in rather general terms 
starting from the properties of the Lie algebra of the generators of the statistical 
density operators as described in ref. [6]. Therefore the HFB approximation is 
only a special case of the recipe described below. 
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2 Determination of the sign of the trace of the sta- 
tistical density operator. 



2.1 Symbols, definitions and basic properties. 

As mentioned in the introduction we shall keep the discussion as general as 
possible. Let Ng be the total number of the single particle states (that is neutrons 
plus protons). Let us consider an arbitrary antisymmetric complex 2Ns x 2Ns 
matrix A and let us define the row vector 7^ = (a, a^), the collection of all an- 
nihilation and of all creation operators. In order to use consistently the matrix 
notations let us denote the column vector 7c = col{a,a^). A general statistical 
density operator (SDO for short) is written as 

W = exp(^7rA7,) (1) 

No other limitations are imposed on this operator, except for the antisymmetry of 
the matrix A. Also let us define the 2Ns x 2Ns matrix 

(2) 

\l 0) 

and the vectors 

7r = IrO- 7c = (^Ic (3) 

Eq. (1) can be recast as 

iy = exp(^7,S7c), S = aA (4) 

This is the form of the SDO we shall use in the following. In the case of the 
HFB approximation S is hermitian. To any operator of the form of eq.(4) one can 



associate the matrix (without the caret) 

W^e^ (5) 

Following ref. [6], the exponents of operators as in eq.(4) form a Lie algebra, and 

therefore the product of any two exponential operators of this form is an operator 
of the same form, moreover the product preserves the association of eq.(5) that is, 
if 

W1W2 = W (6) 

also 

WiW2^W (7) 

For example, Wi is the rotation operator in terms of the Euler angles. In general, 
an arbitrary SDO is constructed as a product of several operators of the class of 
eq.(4), and for each one of them we know unambiguously the matrix S of eq.(4). 
Eq. (6) tells us that the matrix S corresponding to the product W exists, but we 
cannot reconstruct this matrix from eq.(7) because of the 27ri ambiguity of the 
logarithm of the eigenvalues of W. This is the source of the sign ambiguity is the 
evaluation of the traces. From ref. [6], the operators of the class (4) transform the 
vectors 7 in the following way 

where W denote the transpose of the matrix W. Moreover the following relation 
holds (cf. ref. [6]) 

aWa = W'^ (9) 

It ensures that the transformed operators in eq. (8) obey the anticommutation 
relations. 
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2.2 The trace in the case the matrix S is known. 



If the matrix S in eq.(5) is known, the grand canonical trace of W can easily be 
evaluated without ambiguities. Let us prove it in the most general case, assuming 
S is known. Consider a SDO T, then 

IV(W) = lY(t-^Wt) = TV(exp(^f-^7,S7ct)) = 'I¥(exp(^7,T-^ST7,)) 

(10) 

where eq.(8) has been used. We shall prove that the matrix T that diagonalizes S 
satisfies eq.(9). 

It is easy to see that the eigenvalues of S come in opposite pairs. In fact, the eigen- 
value problem for S = a A written in the form (T is the matrix of the eigenvectors 
and the A's are the eigenvalues written in block form for convenience) 



ST^T 



fx o\ 



(11) 



can be rewritten as 



Aa{aTa) = {aTa){a 



A 



a) = {aTa) 



\ 



A' 
A 



A' 

but Aa — —S hence A = —A'. The trace can now be trivially evaluated and the 
result is 

Tr(W) =e-2Si^^'n(l + ^^') (12) 

i 

provided one can show that T is an element of the class of eq.(4), that is, provided 
its associated matrix T satisfies eq.(9) (which guarantees the legitimacy of the 
chain of steps in eq.(lO)). In order to see this, consider the eigenvalue problem for 
W written as 



WT = T 







e 



-A 



(13) 



{Wy^aTa = aTa 



(14) 



W{afa)-^ = {afa)-^ 



\ 



(15) 



by applying at the left and at the right the matrix a ( a^ — 1) and taking into 
account eq.(9) one has 

^ 



Taking the inverse and the transpose of the above we obtain 

^ 
e"^ 

Therefore Tj^ and {aTa)~^ coincide apart a normalization constant 7^ 0. If / 
is the diagonal matrix of elements then 

/ = {afa)T (16) 

Evaluating the above product one can show that / has the doublet structure / = 
diag{fi, /2, .../i, /2, ■■)• Therefore aj'^l'^a — /"^^^ and from this it follows that 
if T does not satisfy eq. (9) the matrix Tj ^/f does. 

Eq. (12) for the trace does not have any sign ambiguity since we had access 
to the eigenvalues of S. If we do not have access to the matrix S, from eq. (12) 
taking the square, we have 

Tr(W)2 = det(l + W) (17) 

since we always have access to eigenvalues of W . This is where the sign ambigu- 
ity comes from. If we work with W obtained from eq.(7), we never have access 
to the matrix S although we know it exists. 

2.3 The contribution of the vacuum. 

Consider now the fugacity dependent trace 

Zg,{z) = Tr(e"^W) (18) 



where z — e'^ and N is the particle number operator. For z = we isolate the 
contribution of the vacuum. Let us define the operator of the class (4) 



and its associated matrix 



a 
-a 

\ 



\ 



J 



7c 



(19) 



(20) 



z 
l/z ) 

Let us set W{z) — J\fW. This operator has an associated matrix W{z) — NW 
Explicitly 

W = 



^ W^ll ^ 



21 



Then 



Zgc{zf = z^^ det(l + W{z)) = det 
This expression can be recast as 



1 

z 



Aet{l +NW) (22) 



Zgc{zY = det(5^ + zSp) = det(5„) det(l + zS'^Sp) 



where 





/ 


1 








W,2 


\ 








'] 








I 


W21 


W22 ) 




[ 


1 


/ 



(23) 



(24) 



Since Zgc{z) is a polynomial in z the eigenvalues of S~^Sp must come in degen- 
erate pairs {nk: A*fe), k — 1, , , A^j. This argument is the similar to the one used in 
ref. [7] .Therefore 



Zg,{z) = det{W22Y^^l[il + ziik) 

k=l 



(25) 
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It follows that, if we know the matrix S (the log of W) then, using eq.(12). 

Where we have set z — 1. There are no sign ambiguities in n£i(l + ^A*fc) since 
its sign must be a continuous function of z. The only ambiguity is the det^*^^ (1^22) 
and it is removed by eq.(12). Eq.(26) is the basic equation that allows us to remove 
the sign ambiguity also in the general case when we do not know the matrix S and 
its eigenvalues unambiguously, as shown in the next subsection. Before leaving 
this section let us note that for z = we obtain the vacuum contribution to the 
Grand Canoncal trace of W 

< Q\W\(} >= det(iy22)^/^ (27) 

where |0 > is the particle vacuum. 

2.4 The trace in the general case. 

Consider the SDO of the type 

W = W^^W^W (28) 

where W^''^"-^ are of the type of eq.(4), and let us assume we know explicitly the 
matrices -S"^^) and -S^"^. As previously mentioned, we do not know unambiguously 
the matrix S associated with W, although we know the matrix W associated with 

W, since W = W''^^W^"-\ In what follows we shall need the following matrices 

= Wi-^wi\ C(«) = Wi^wif-' (29) 
Using eq.(25) (with 2; = 1) we have 

= {det[W^'^W^%2y/' ft (1 + /^^) (30) 

k=l 
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The square root of the determinant can be evaluated in the following way. Con- 
sider the vacuum expectation value < 0|iy('')e"'^W*^'^)|0 >. Then (cf. eq.(27)) 

Direct evaluation of the matrix product gives 

or 

< 0|iy(^)e"^W(^)|0 >= (detwgV/'(detwgV/2(det(l +z2d(^)C(^))V2 

(33) 

Since this vacuum contribution must be a polynomial in z, the eigenvalues of 
]j{b)(j{a) must come in degenerate pairs {vk,i^k)- Therefore, considering only one 
eigenvalue for each degenerate pair, 

< 0|W^(^)e"^W(^)|0 >= (det wgV/'(det wg)i/2 -Q^l ^ ^2^^^ ^34^, 

k=l 

Finally setting z — 1, we obtain for the grand-canonical trace 

Z,, = det(w4'V/' det(iyi"V/' 11(1 + Uk) 11(1 + i^u) (35) 

k k 

The only sign ambiguity in eq.(35) comes from the contributions of the two square 
roots. From eq.(34) (for z = 0), one can see that each square root is again the 
vacuum contribution from W^'^^ and W^"'^ but we know already as to remove this 
ambiguity using eq. (26) for each W'^^'> , 1^*^"^, since we know the matrices 5'^'') and 
5"*^"^ These considerations can easily be extended to a product of several SDO's. 

As a final remark, it is possible (using eq. (9)) to prove that the matrices D^''^ 
and C*^") are antisymmetric and therefore these arguments amount to a quantum 
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mechanical proof of the statement that the product of two antisymmetric matrices 

has eigenvalues in degenerate pairs (even A^^) or the odd one is zero (odd Ng). 

2.5 A numerical test. 

Essentially our method to fix the sign of the square root in the general case 
is based on an analytical continuity argument, supplemented by the fact that we 
know the contribution of the vacua of the factors W^"-\W^^\ .., because of eqs.(12) 
and (26). We performed extensive numerical tests of eq.(35), by considering an 
ensemble of 12 x 12 antisymmetric random matrices with matrix elements uni- 
formly distributed in the interval —a and a (a ~ 3). This random set generates the 
matrices A^^'> and A'^"-'' and from eq.(4) the matrices S^^'> and S^"'\ To test the above 
method we consider the matrices hS^^^ and hS^"'^ with h varying from to 1 in suf- 
ficiently small steps, so that a numerical continuity argument can be tested. We 
also have used the conventional way of evaluating the grand canonical trace, by 
taking the square of eq.(35) and then numerically evaluating the square root and 
then following the phase of this square root as h is varied from to 1. More pre- 
cisely, we evaluate the phase of the grand canonical trace and, if the phase changes 
between consecutive values of 6 by an amount larger than some value 5(t)max we 
change the phase by tt. In many cases the numerical continuity argument agrees 
with eq.(35) but in some case we found a sign disagreement. It is instructive to 
analyze these latter cases. In fig.(l) we show the phase of eq.(35) as a function of 
h for one of these instances, and in fig. (2) we show the phase for the same trace 
evaluated with the numerical continuity argument with two different step sizes. 
For a relatively small number of steps {N steps — 200), the phase shows a discon- 
tinuity, while for a much larger number of steps (Ngteps — 800), the continuity of 
the phase as a function of b is restored. The maximum phase change from one 
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Figure 1: The phase of eq.(35) as a function of b 

step to the next has been kept fixed for both cases to S(f)max — 0.5. From fig.(2) 
we see that the source of the discontinuity is the rather large change in the phase 
for small variations of b, which is associated with a vertical slope of the phase. 
In this instance even 400 steps fail to reproduce the continuity of the phase. This 
rather surprising result, was obtained because of the disagreement with eq.(35), 
otherwise it would have gone undetected. If we have to check the continuity of the 
phase for hundreds of cases (as it is in the case of the angular momentum projec- 
tion), we can hardly check every single instance to insure the proper phase. If we 
consider smaller values of 50maa; we may need several thousands steps to restore 
continuity. This pathological behavior of the numerical continuity argument was 
found in presence of a vertical slope . This example does show the limitations of 
the numerical continuity method. 
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3 Conclusions. 



In this work we have shown that the sign ambiguity in the evaluation of the trace 
of statistical density operator written as a product of elementary statistical density 
operators (this is usually the case in physical applications) can be effectively re- 
moved by considering all factors separately, evaluating the vacuum contributions 
using eq.(12), then reconstructing the vacuum contribution of the full statistical 
density operator without taking any square root of determinant using eqs. (35). 
Quite surprisingly, a numerical continuity argument, computationally more in- 
volved, can fail to reproduce the proper sign in some numerical test cases, unless 
we use a large number of steps. The recipe presented in this work opens the pos- 
sibility to perform calculations of grand canonical partition functions within the 
HFB formalism with projectors to good quantum numbers. 
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